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ABSTRACT 

A major recent development in extragalactic astronomy has been the discovery of a 
population of galaxies that are luminous at submillimetre wavelengths. Estimates of 
their spectral energy distributions suggest that these galaxies are the high-redshift 
analogues of the ultraluminous infrared galaxies (ULIRGs) observed locally. Here we 
investigate the implications for the galaxy formation problem if this is indeed the case. 
We identify plausible (but non-unique) redshift-dependent galaxy luminosity functions 
that are consistent with both source counts at 2800, 850, 450 and 175 ^m, and far- 
infrared background radiation intensities at 850, 240 and 140 /im. In all our models, 
most of the submillimetre-luminous sources are distant galaxies with high bolometric 
luminosities > 1O 12 L . As for many local ULIRGs, it is not possible to determine 
whether these luminous galaxies are powered by starbursts, like the local galaxy Arp 
220, or by active galactic nuclei (AGN), like the local galaxy Markarian 231. We 
investigate both possibilities. If the submillimetre-luminous galaxies are all starbursts, 
then we predict that the fraction of the cosmic star-formation rate in these objects 
is large, but does not necessarily dominate the star-formation rate of the Universe - 
the Madau plot - at any redshift. Only a few per cent by mass of the present-epoch 
spheroidal stellar population would have been formed in such a population of star- 
forming galaxies, consistent with the constraints on the number of galaxies with old 
stellar populations in the field at low and intermediate redshifts derived from A"-band 
surveys. If the submillimetre-luminous galaxies are all powered by AGN, then the 
comoving density of supermassive black holes onto which material is accreting at high 
redshift probably equals no more than a few per cent of the local density of massive 
dark objects. 

Key words: galaxies: formation - galaxies: evolution - infrared: galaxies - diffuse 
radiation - quasars: general - cosmology: observations 



1 INTRODUCTION 

Recent observations with the Submillimetre Common-User 
Bolometer Array (SCUBA; Holland et al. 1999) on the 
James Clerk Maxwell Telescope have highlighted the pres- 
ence of a number of submillimetre-luminous galaxies (Smail, 
Ivison & Blain 1997; Holland et al. 1998; Barger et al. 1998a; 
Hughes et al. 1998; Eales et al. 1998). To date about forty 
sources have been found. These measurements confirm ear- 
lier suggestions (Blain & Longair 1993; Eales & Edmunds 
1996, 1997) that submillimetre-wave observations will pro- 
vide an important probe of cosmology. 

The optical counterparts of these SCUBA sources are 
faint, most with J > 20 (Smail et al. 1998). They are 
presumably not low-redshift sources: the best studied is 



SMMJ02399-0136 (Ivison et al. 1998), a dust-enshrouded 
AGN at z = 2.8, perhaps similar to IRAS F10214+4724 
(Rowan-Robinson et al. 1991; Lacy et al. 1998). The recent 
detection of molecular gas in SMM J02399-0136 (Frayer et 
al. 1998) also suggests that these galaxies are similar. 

What is the nature of these SCUBA-selected sources? 
In particular are they similar to the ultraluminous infrared 
galaxies (ULIRGs) observed locally (e.g. Sanders & Mirabel 
1996)? They have the same submillimetre-wave flux den- 
sities as would have local ULIRGs seen at high redshift 
(Barger et al. 1998a). They also have similar optical colours 
(Smail et al. 1998) to two of the three local ULIRGs stud- 
ied at ultraviolet wavelengths by Trentham, Kormendy & 
Sanders (1998), were those galaxies to be seen at high red- 
shift. If this interpretation is correct, then we might hope to 
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be able to use the wealth of observational information avail- 
able for local ULIRGs to help to understand the properties 
of the SCUBA sources and their relevance to cosmology. 

Whether the SCUBA sources are dust-enshrouded 
AGNs, like Markarian 231, or dust-enshrouded starbursts, 
like Arp 220, is an important question. This is difficult to 
determine for even local ULIRGs, because there are up to 
several hundred magnitudes of extinction along our lines of 
sight to the galaxy cores at optical wavelengths. For local 
ULIRGs, simple conclusions can be drawn based on the 
general form of the spectral energy distributions (SEDs) 
(e.g. Sanders et al. 1988a); however, this information is not 
available for the SCUBA sources. Recently, new more de- 
tailed methods based on mid-infrared spectroscopic diag- 
nostics (Lutz et al. 1996, Genzel et al. 1998) have been used 
to resolve this question for the most local ULIRGs. 

We construct redshift-dependent luminosity functions 
that are consistent with all the counts and backgrounds 
from a number of recent surveys at far-infrared and millime- 
tre/submillimetre wavelengths. We use non-evolving Gaus- 
sian luminosity functions over specified redshift ranges. 
These are the simplest luminosity functions that we could 
adopt, without the results being dependent on the prop- 
erties of low-luminosity galaxies that are not probed by 
the SCUBA observations. Contributions from these low- 
luminosity galaxies are important when computing back- 
grounds, although they contribute insignificantly to the 
counts. The parameterization used to compare models with 
observations reflects this fact. Imposing all the constraints 
simultaneously does limit the possible form of the luminosity 
function and we identify three plausible models. 

Given the redshift-dependent luminosity functions of 
these models, we investigate the properties to the individ- 
ual sources using observations of their local ultraluminous 
counterparts, and derive the cosmological implications of the 
SCUBA galaxies under both the starburst and AGN inter- 
pretation. 

We first investigate the possibility that the SCUBA 
sources are high-redshift star-forming galaxies. The star- 
burst models of Leitherer & Heckman (1995) are used to 
convert far-infrared luminosities to star-formation rates; a 
transformation that is uncertain by more than an order 
of magnitude. We then place the SCUBA galaxies on the 
Madau plot, which relates the cosmic star formation rate to 
the cosmic epoch (Madau et al. 1996; Madau, Delia Valle & 
Panagia 1998), and predict the total density of local stars 
produced in such objects. Most of the star formation in 
the SCUBA sources is observed through huge amounts of 
internal extinction, and so will not be included in global 
star-formation rates that are computed using the optically 
selected samples that are normally used to construct the 
Madau plot. Indeed Hughes et al. (1998) show that deriv- 
ing a star-formation rate from a rest-frame ultraviolet flux 
results in a value that is more than an order of magnitude 
too low. The SCUBA sources are thus not accounted for in 
the existing samples that are used to construct the Madau 
plot. 

We then investigate the possibility that the SCUBA 
sources are high-redshift dust-enshrouded AGN, as dis- 
cussed by Haehnelt, Natarajan & Rees (1998) and Almaini, 
Lawrence & Boyle (1999). The three distant submillimetre- 
luminous galaxies that have been studied in detail - 



APM 08279+5255 (Irwin et al. 1998; Lewis et al. 1998; 
Downes et al. submitted), SMM J02399— 0136 and IRAS 
F10214+4724 - all contain powerful AGN. If all the sources 
derive their bolometric luminosity from AGN that heat their 
dust shrouds radiatively by accretion onto a massive black 
hole, then we can use our redshift-dependent luminosity 
functions to compute the comoving integrated mass density 
of these black holes. 

We find plausible (although not unique) scenarios that 
are consistent with all the observations and present their 
cosmological implications. Finally, we highlight future work 
which may help to distinguish between the various scenar- 
ios. Throughout this paper we assume an Einstein-de Sitter 
world model with Ho = 50kms~ Mpc -1 . 



2 FAR-INFRARED /SUBMILLIMETRE-WAVE 
COUNTS AND BACKGROUNDS AND THE 
ULIRG LUMINOSITY FUNCTION 

There have been many recent measurements of the far- 
infrared and submillimetre source counts and backgrounds 
at a number of wavelengths: see Table 1. This has been a 
recent field of substantial activity because of new instrumen- 
tation, both SCUBA and ISO. If we assume that the SED at 
far-infrared and submillimeter wavelengths is known for all 
the sources, we can then use all these measurements in con- 
junction to constrain the bivariate far-infrared luminosity- 
redshift function of these sources. 

We assume that the galaxies have the thermal SEDs 
of warm dust that is heated by an enshrouded starburst or 
AGN, and adopt a simple analytic form for the luminosity 
function between a minimum and maximum redshift. All the 
observational constraints on the model are imposed through 
integral functions, and so do not require unique solutions. 
Nevertheless, it is interesting to see which general classes of 
luminosity-redshift distributions are ruled out and why. 

To ensure consistent definitions are used when matching 
models to observations, first we present all the details of 
our computations. Much of this will be familiar to many 
readers, but different authors define parameters in slightly 
different ways. Secondly, we point out some generic features 
of the comparison between our models and observations, in 
particular the results of requiring the models to be consistent 
with source counts at one wavelength and with the infrared 
background at another simultaneously. Finally, we isolate 
some plausible luminosity functions that are consistent with 
all the observations for further investigation. 

We define a bivariate 60- /im luminosity-redshift 
function 4> z (Leo), with units Mpc _3 Lg 1 , such that 
<?^z(i6o) dZ/6o dz is the total number density of galaxies with 
60-pim luminosity between Leo and L6o + dL6o with redshifts 
between z and z + dz. Some familiar analytic examples of 
this function are a Gaussian in log 10 Leo with no luminosity 
or density evolution, 



Z (L 60 ) = C(l + «) 3 exp 



1,2/' Leo 
- 2^2 logic ( — 



1 



Leo lnlO 



(1) 



or a Saunders et al. (1990) function, which is a power-law in 
L 6 o with index a at the faint end and a Gaussian in log 10 Leo 
at the bright end, and allows for density and luminosity 
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Table 1. Far-infrared and submillimetre-wave surveys. Lagache et al. (1998) and Kawara et al. (1997, 1998) obtain a different calibration 
for the ISO 175-^tm counts. A discussion of calibration is given in Lagache et al. (1998). The detection limit in the 850-/im survey by 
Hughes et al. (1998) is about 2 mjy: the count at 1 mjy is derived from a source confusion analysis. See Blain et al. (1999b) and Smail et 
al. (1999) for a direct sub-mjy 850-/im count. For complementary far-infrarcd/submillimetre-wave background measurements the reader 
is referred to Puget et al. (1996), Guiderdoni ct al. (1997), Dwck et al. (1998) and Fixsen et al. (1998). 



Wavclc ngt h 


Counts / 
deg~ 2 


Flux Limit / 
mjy 


Background / 
nW m sr _1 


Telescope /Instrument 


Reference 


2.8 mm 


< 162 


3.55 




BIMA 


Wilner & Wright (1997) 


850 Mm 


2500 ± 1400 
1100 ± 600 
800™ 

7000 ± 3000 
1800 ± 600 


4 
8 
3 
1 

2.8 


— 


JCMT/SCUBA 


Smail ct al. (1997) 
Holland ct al. (1998) 
Barger et al. (1998a) 
Hughes et al. (1998) 

Ealcs ct al. (1998) 


850 /xm 






0.55 ±0.15 


COBE/FIRAS 


Fixsen et al. (1998) 


450 /im 


< 1000 

< 360 


80 
75 




JCMT/SCUBA 


Smail ct al. (1997) 
Barger et al. (1998a) 


240 /xm 






17 ±4 
14 ±4 


COBE/BIRBE + IRAS/ISSA 


Schlegel et al. (1998) 
Hauser et al. (1998) 


175 Mm 


41 ±6 
98 ±15 


150 
100 




/SO/ISOPHOT 


Kawara ct al. (1998) 
Lagache et al. (1998) 


140 Mm 






32 ± 13 
24 ±12 


COBE/BIRBE + IRAS/ISSA 


Schlegel et al. (1998) 
Hauser ct al. (1998) 



evolution, 



in units of sr 1 (or deg 2 ), where, 



^60 



exp 



2a 2 



lo Sio 



60 00 
1 + 



Leo In 10 



(2) 



Implicit in both equations (1) and (2) above are a normal- 
ization constant C (units Mpc~ 3 , which depends on z if 
there is density evolution), a characteristic luminosity Lg 
(which depends on z if there is luminosity evolution), and a 
Gaussian width a. 

Once we have specified 4> z (Leo) for some population of 
galaxies, we can compute their contribution to the cosmic 
infrared background at some frequency v: 



Leo) 



vl _ i r r l(z, 

^ Jo Jo 4nd L (z)2 v 
in units of Wm~ 2 sr -1 , where, 



dV 

, (L 6 o) dLeo -3— dz, 
dz 



l v (z, Leo) = 



veo 



Leo 



exp(^sa* 



kT 



exp [ 



kT J 



(3) 



(4) 



and veo is the frequency corresponding to a wavelength 
of 60 lira. k v is the emissivity function of dust. We follow 
Hughes (1996) and Blain et al. (1999a) in assuming a power 
law k v ~ f 1 ' 5 , and so the dust emission spectrum at long 
wavelengths is a Raleigh-Jeans power-law with spectral in- 
dex 3.5. We can also compute the number density of sources 
with flux density above some threshold Sum, which is mea- 
sured in units of Wm -2 Hz -1 (or Jy), 



n(Si im ) = — 



ilim[Sn m ] 



4>z(Leo) dLeo ^p- dz, (5) 
dz 



£lim[Sl im ] =4ndl J (z) 2 Slim ^60 



veo 



1/(1 + z) 

Thvt 

exp[— \ 



l±£l " 



The quantities, 



di,(z) 

and, 

dV 

dz = 



2c 
~H~o 



1 - 



k v(1+z) exp(^a)-l 



(l + z), 



is. G|) 3 (i + o- 



1 - 



(6) 



(7) 



(8) 



are the luminosity distance and the redshift-derivative of the 
volume clement respectively (assuming fio = !)■ 

A few features of the comparison between observation 
and theory are immediately apparent. Source counts con- 
strain the integral of 4> z (Leo) over redshift and luminos- 
ity above some redshift- dependent limit. However, measure- 
ments of the background light constrain the integral of the 
product of luminosity Leo and 4> z (Leo) over redshift and 
all luminosities. For a Saunders luminosity function with 
pure density evolution, wc compare in Fig. 1 the ratio of 
the predicted 850-/im source count and 240-/im background 
as compared with the observed values listed in Table 1 as 
a function of the characteristic luminosity L* = Lg . The 
minimum is produced by the different constraints on the 
integral of <f> z (Leo) imposed by the source counts and the 
background measurement. 

Very low luminosity galaxies make a negligible contri- 
bution to the source counts, but a substantial one to the 
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Figure 1. The normalization obtained by fitting a Saunders 
function (equation 2) with density evolution parameterized by 
7 to the COBE 240-^tm background (Schlegel et al. 1998), rela- 
tive to the normalization obtained by fitting a similar function 
to the SCUBA 850-/mi counts of Smail et al. (1997). A value of 
a = 0.724, derived by Saunders et al. (1990), is assumed. A dust 
cmissivity index of 1.5 is assumed when converting luminosities 
measured at different wavelengths to restframe 60-/j.m luminosi- 
ties. The four line styles represent the following parameter values: 
solid (T = 70K, 7 = 0), short-dashed (T = 40 K, 7 = 0), long- 
dashed (T = 70 K, 7 = 6), and dot-dashed (T = 40 K, 7 = 6). 
Locally, L* m = 1.1 X 10 9 L . 



infrared background if L* is low in the Saunders function. 
Hence, a low normalization is required to explain the back- 
ground data, but a high one is required to explain the source 
count data, and so the ratio plotted in Fig. 1 is very high 
at low L». Increasing L* increases the fraction of high- 
luminosity galaxies, and so lowers the ratio of the normal- 
izations required in the figure. At very high L„, the source 
counts are produced by a very small number of extremely lu- 
minous sources, which cause the background to be extremely 
high based on their luminosities alone. Hence, the ratio of 
the normalization begins to increase again, producing the 
minimum in Fig. 1. If the flux threshold (lower limit) in the 
luminosity integral in equation (5) is reduced, then the turn- 
up of the ratio at high luminosities is less marked, and dis- 
appears as the lower limit tends to zero. 

It is intriguing that this minimum in the ratio occurs 
at a very high characteristic luminosity L* = 10 12 Lq. Lo- 
cally L, ~ 10 9 Lq , and so if the luminosity function of the 
SCUBA sources has a Saunders form, huge amounts of lumi- 
nosity evolution are required to match the data. This was es- 
sentially one of the main conclusions of Blain et al. (1999a). 

Higher temperature sources systematically produce 
greater background fluxes and counts at shorter wave- 
lengths, because their dust spectra peak at shorter wave- 
lengths. Hence, lower normalizations are required in Fig. 1 
in order to explain the data for 40-K sources as compared 



with 70-K sources. The 40-K curve is consistent with a ra- 
tio of one over a range of luminosities, suggesting that if a 
Saunders luminosity function describes the SCUBA sources, 
then they must have temperatures of about 40 K. This was 
another conclusion of Blain et al. (1999a): see their Figure 
4 - in their favoured scenarios most of the SCUBA sources 
are at z > 2 and have luminosities ~ 10 12 L Q and dust tem- 
peratures ~ 40 K. 

Finally, note that the results in Fig. 1 are only weakly 
dependent on the form of density evolution parameterized 
by 7. 



2.1 Comparison with the 850-pim source counts 

We now compute the normalizations C required to fit the 
850-pm source count data for a number of Gaussian lumi- 
nosity functions (equation 1) with no density or luminosity 
evolution. C is constant between redshifts z\ and Z2 and 
is zero elsewhere. This is the simplest possible form of the 
luminosity function; at present the observations do not jus- 
tify a more rigorous treatment. This parametrization is mo- 
tivated in part by the fact that the local 60-/j,m luminos- 
ity function (Saunders et al. 1990) is approximately Gaus- 
sian at the bright end. The local elliptical galaxy luminos- 
ity function (Binggeli, Sandage & Tammann 1988, Ferguson 
& Sandage 1991) is also Gaussian, which is relevant if the 
SCUBA sources evolve into elliptical galaxies: see Section 
3. One important feature of the Gaussian function is that 
it is only valid for high-luminosity galaxies. Large numbers 
of galaxies with low far-infrared luminosities will also exist 
between z\ and 22 . Although these will not contribute to the 
source counts, they may contribute significantly to the in- 
frared backgrounds if they are very numerous. Hence, when 
we match our luminosity functions to counts (equation 5), 
we require the model prediction to equal the observed values, 
but when we match our luminosity functions to background 
fluxes (equation 3), we require the model prediction to be 
less than or equal to the observed values. Our approach is 
fundamentally different from that of Blain et al. (1999a). 
We are not evolving the local 60-/xm luminosity function 
to high redshift. Indeed, in Section 3, we argue that the 
systems which have the highest far-infrared luminosities for 
redshifts z\ < z < Z2 probably evolve into elliptical galax- 
ies, which are not the most luminous systems at far-infrared 
wavelengths in the local Universe. 

In Figs 2 to 7 we present the required normalization con- 
stants for a Gaussian luminosity function to fit the 850-/im 
data of Smail et al. (1997) for various values of a and L* 
given a constant dust temperature T. The uncertainties are 
large (see Table 1), and so the range of acceptable normaliza- 
tions is large, with a peak-to-peak spread of about a factor 
of 4. Given these large uncertainties, all the 850-/im counts 
by different authors listed in Table 1 are fully consistent. 
Some general features of Figs 2 to 7 arc worth highlighting: 

(i) (j = (a (^-function) never appears to fit the data well; 

(ii) lower normalizations are required for higher charac- 
teristic luminosities; 

(iii) the absolute values of the normalizations depend pri- 
marily on where the dust spectrum peak is shifted in wave- 
length space. It is shifted closer to 850 /u,m for a source at 
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Figure 2. Combinations of the normalization constant C and the width a for Gaussian luminosity functions with no density 
evolution that are consistent with the 850- (im counts. The units of C are Mpc - 3 . All galaxies arc assumed to have a temperature 
T = 70 K. The comoving galaxy density C is assumed to be constant between rcdshifts z\ = 3 and 22 = 5 and zero outside this 
range, (that is the upper and lower limits to the redshift integral in equation (5) are z\ and 22)- The solid lines show fits to the 
mean value quoted by Smail et al. (1997). The dottcd-dashed lines are fitted to the +lc values. The dashed lines are fitted to the 
— values: see Table 1. The four panels are plotted for different values of L». 



higher redshift at either 40 or 70 K, and so a lower normal- 
ization is required at higher 2; 

(iv) lower temperatures require lower normalizations for 
a given 60-^im luminosity, since 850 /im is closer to the black- 
body peak. 

2.2 Comparison with counts at other wavelengths 

Having isolated various models that are consistent with the 
850-^im counts, we now consider the constraints from mea- 
surements at other wavelengths. We select a number of mod- 
els with various values of a, L», C, and T, each represented 
by a point on the curves in Figs 2 to 7, and thus consistent 



with the 850-pim count. For each model, we consider three 
separate scenarios: scenario "0", in which the 850-/im source 
counts are the mean Smail et al. (1997) values; scenario "+", 
in which they are la larger; and scenario "— ", in which they 
are la smaller. The counts in the "— " scenario are very close 
to the lowest published counts (Barger et al. 1998a). In all, a 
total of 45 models, each of which generates three scenarios, 
are considered. 

In Table 2, we show how well the models defined by 
these parameters account for the measurements at other 
wavelengths. The numbers in that table represent the model 
prediction relative to the observed counts or background 
fluxes listed in Table 1, assuming scenario "0". The frac- 
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tional uncertainties in the numbers in Table 2 are given in 
Table 3. 

Other measurements exist in the literature that we do 
not use while making this comparison. For example, mea- 
surements of 15-/im source counts with ISO have recently 
been made (Rowan- Robinson et al. 1997, Aussel et al. 1998). 
However, these measurements probe far down the Wien tail 
of the dust spectrum, and are of very limited use in con- 
straining the high-redshift luminosity function. Furthermore 
we would need to account for contamination in the 15-pim 
samples by nearby (unobscured) galaxies and AGN. Simi- 
larly the IRAS 60-^im source counts (Saunders et al. 1990) 
are dominated by nearby (z < 0.1) sources. Measurements 
at 1.25 mm offer a promising probe in the future, but to 
date measurements have only been made for IRAS-selected 
galaxies (Franceschini, Andreani & Danese 1998). Unbiased 
samples are not yet available at this wavelength. 



Some of the more notable general trends shown in Table 
2 are: 

(i) models in which the sources are at higher redshift 
tend to produce higher source counts or backgrounds at 
longer wavelengths, because the peak of the dust spectrum 
is shifted to longer wavelengths; 

(ii) models in which the sources are at higher redshift 
tend to produce higher backgrounds for a given source count 
at the same wavelength. The counts probe only part of the 
total range of galaxy luminosities, and there are a larger 
number of low-luminosity sources that contribute to the 
backgrounds but not to the counts; 

(iii) on decreasing the temperature from 70 to 40 K, the 
predicted 2. 8- mm counts increase, but the predicted 175-/im 
counts and 140- and 240-pim backgrounds decrease, because 
the peak in the dust emission spectrum shifts to longer wave- 
lengths; 
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(iv) Changing the temperature has a substantial effect 
on the relative value of the source counts and backgrounds 
at the same wavelength. Decreasing the temperature has a 
greater relative effect on the counts, because the lower limit 
in the integral in equation 5 is decreased: see Fig. 8. 

Seven of the models that we investigate are consis- 
tent with all the observations: see the last column of Ta- 
ble 2 for details. These all have a high characteristic lumi- 
nosity greater than 1O 11 L0. Using a different parametriza- 
tion, Blain et al. (1999a) also found that most of the 
SCUBA sources are distant galaxies with similarly high lu- 
minosities. These seven models are: 12/40/3/5/0. 9"0" (here- 
after Model A); 13/40/3/5/0.9"-" (hereafter Model B); 
11/40/2/4/0.9"-" (hereafter Model C); 12/40/2/4/0.9"-" 
(hereafter Model D); 13/40/2/4/0.5"0" (hereafter Model E); 
13/40/2/4/0.5"-" (hereafter Model F); 13/40/1/3/0.1"-" 
(hereafter Model G). In general, the tables suggest that 



models with T — 70 K overproduce the 175-^im counts and 
the 140- and 240-/im backgrounds, given the normalization 
from the 850-/xm source counts. Models in which the galax- 
ies are nearby, that is where z\ = 1, tend to have similar 
problems. Conversely, models with very low temperatures 
T < 40 K would overproduce the 2.8-mm counts, but such 
low-temperature sources are ruled out by the consistency 
arguments of Blain et al. (1999a — see Section 4 of that 
paper). 



2.3 Comparison with the FIRAS background 

Fixsen et al. (1998) derived the cosmic submillimetre back- 
ground from COBE FIRAS data. Although the background 
at these wavelengths is dominated by emission from the 
Galaxy and the cosmic microwave background, they used 
three independent techniques to subtract out these signals. 
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Figure 5. As Figure 2, but for T = 40K. 



The resulting residual 850-/im extragalactic background ra- 
diation intensity ^/„=0.55 ± 0.15 nWm~ 2 sr _1 . 

We can predict 850-/im backgrounds in our models us- 
ing equation (3). The results are presented in Table 4. As 
in the previous section, we require that our models do not 
overproduce the background. Model G achieves this con- 
vincingly, and models B and F are consistent within 2a. We 
select these three models for further study. In general, mod- 
els with «i = 3 overpredict the 850-pm background if the 
240-^im background is predicted correctly, because the ob- 
served spectra at a temperature of 40 K are shifted to too 
long a wavelength. The surviving models have z\ = 2 (B and 
F) or 2i = 1 (G). All these models also have L* = 10 13 L Q , 
meaning that most of the far-infrared luminosity from this 
population comes from galaxies with 60-^im luminosities in 
excess of 10 12 Lq. As discussed in Section 2, models with 
lower characteristic luminosities tend to overpredict the 850- 
fim background if they match the 850-pm counts. 



Note that in the three surviving models (B, F and G) 
the normalization of the 850-^im counts is consistent with 
all the observations. In the best-fitting model (G) , the mean 
redshift of the sources z B ~ 2, as argued by Lilly et al. (1998) . 
One half (F,G) or more (B) of the 240-/im background, but 
considerably less of the 140-^im background, comes from the 
high- luminosity galaxies described by our Gaussian luminos- 
ity function. In models B and F, about half of the 175-/im 
sources are the same sources that contribute to the 850-/im 
source counts; in model G, the two populations are practi- 
cally the same. In Model B, the predicted 2.8-mm count is 
close to the observed limit (Wilner & Wright 1997); for the 
other models the predicted count is much smaller. The 450- 
/jm source count limit (Barger et al. 1998a) is larger than 
the predicted values by a factor of about 3 for Model B, and 
by a much larger factor for the other models. 

Three models appear to fit all the data. Even within our 
simple parametrization we do not find a unique best-fitting 
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Figure 6. As Figure 3, but for T = 40K. 



model. Furthermore, there are presumably many other (non- 
Gaussian) models that fit all the data, for example the mod- 
els of Blain et al. (1999a). Nevertheless most of our models 
seem to be ruled out when we consider all the observations 
in conjunction. That a relatively narrow region of param- 
eter space is consistent with observation (L* ~ 10 1 Lq, 
T ~ 40 K, z s ~ 2) is encouraging and suggests that it will 
be productive to investigate these models in further detail. 



3 PROPERTIES OF THE SCUBA SOURCES 

We have isolated three models of the redshift-dependent lu- 
minosity function of ULIRGs that are consistent with obser- 
vation, and now investigate the cosmological implications. 
In Section 3.1 we investigate the possibility that these are 
dusty star-forming galaxies. We use the star-formation mod- 
els of Leitherer & Heckman (1995) to place the galaxies on 



the Madau Plot, and then examine the consequences of this 
interpretation in the context of the wider galaxy formation 
picture. In Section 3.2 we investigate the alternative possi- 
bility that the sources are dust-enshrouded AGNs. 



3.1 The SCUBA sources as star-forming galaxies 

The majority of local ULIRGs appear to be star-forming 
galaxies and not dust-enshrouded AGN (Sanders et al. 
1988a, Genzel et al. 1988). Furthermore, the temperatures 
of the SCUBA sources inferred in Section 2 (T ~ 40 K) arc 
close to, although systematically slightly cooler than, those 
of local star-forming ULIRGs. Locally, dust-enshrouded 
AGN like Mrk 231 tend to have higher dust temperatures 
than dust-enshrouded starbursts like Arp 220 (Sanders et at. 
1988a); however, note that some dusty high- redshift quasars 
appear to be fairly cold: see Benford et al. (1998). Hence, 
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it seems reasonable to investigate the possibility that the 
SCUBA sources are high-redshift star-forming galaxies. 



3.1.1 Leitherer-Heckman models 

In order to compute the cosmological star-formation rate as- 
sociated with the SCUBA sources at high redshift, we need a 
recipe to convert far-infrared luminosities to star-formation 
rates for individual sources. 

The transformation between these two quantities is not 
straightforward to determine observationally even for local 
ULIRGs because of the high internal extinction along our 
lines to the galaxy centers. For example, even the near- 
infrared 2.2-/im Br7 line strength-far-infrared luminosity 
correlation (Goldader et al. 1997a, b), which is usually very 
reliable, breaks down in ULIRGs at very high luminosities, 
presumably due to internal extinction. We therefore need 



to compute this transformation using models (Leitherer & 
Heckman 1995). 

In Table 5 we list a number of star-formation mod- 
els that we investigate in further detail. Note that we use 
capital letters to represent galaxy formation models (Sec- 
tion 2) and roman numerals to represent the different star- 
formation models. There is insufficient information to de- 
termine accurately the appropriate set of input parameters 
to the star-formation models for the SCUBA sources, and 
even for local ULIRGs in many cases. The more important 
parameters are: 

(i) the mode of star formation. For example, the star for- 
mation may occur continuously, as an instantaneous burst, 
or with a more complex time dependence, for example 
an exponential decay. We investigate both the instanta- 
neous (Models I and II) and continuous (models III through 
X) cases. Observations of the core of Arp 220 (Mouri & 
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Table 2. Comparison between observations and models. All the numbers represent the predicted source counts or background flux for 
model "0" that are appropriate to the luminosity function defined by the parameters L„, T, a, zi, and Z2, relative to the observed 
counts and background intensities listed in Table 1. The model parameters are listed in condensed form in the model name, which takes 
the general form logio(Lg )/T/ zi/ Z2/0-. The figure panel in which the constraints for each model arc imposed are listed in the second 
column. A value less than 10~ 10 is listed as ~ in the Table. The errors in these quantities are listed in Table 3. All the models have 
been selected a priori to give the correct 850-^tm counts. The comments refer to the selection or rejection of models based on the listed 
values. The range of acceptable values and a description of the numerical code values are given in Table 3. The 140-/im and 240-^tm 
background values of Schlegel et al. (1998), the 175-^tm counts of Lagachc ct al. (1998), and the 450-^tm counts of Smail et al. (1997) 
are fitted. 
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Table 3. Fractional uncertainties in the background and source counts values used to accept or reject the models listed in Tabic 2. The 
acceptable ranges of values are used to make the selection, and the relevant cases are described by numerical codes in the last column 
of Table 2. A code is shown if the model is rejected on the grounds of: (1) overproduction of the 2.8-mm count; (2) overprcdiction of 
the 450-/im count; (3) overproduction of the 175-/^m count; (4) underprediction of the 175-/^m count; (5) overprediction of the COBE 
240-^im background; and (6) overprediction of the COBE 140-/im background. Although a model that undcrprcdicts the 175-^m count 
is not formally excluded, such a model requires that the 175-/im counts and 850-/im counts come from entirely different populations 
of galaxies. This is not impossible, but we regard it as unlikely; we choose a factor of ten discrepancy as the cutoff between a model 
being acceptable and unacceptable. In the last column of Table 2, superscripts 0, + or — indicate that the exclusion applies only in that 
scenario, corresponding to the values listed below. The absence of a superscript implies that the model is rejected in all three scenarios. 



Model 2.8-mm 450-/im 175-/im 240-^tm 140-//m 

counts counts counts background background 



Fractional Error in ratio 






0.15 


0.24 


0.41 


Range of Acceptability ("0" models) 
Range of Acceptability ( "— " models) 
Range of Acceptability ( "+" models) 


< 1.0 

< 2.4 

< 0.63 


< 1.0 

< 2.4 

< 0.63 


[0.85,1.2] 
[2.0,2.8] 
[0.54,0.73] 


< 1.2 

< 3.0 

< 0.78 


< 1.4 

< 3.4 

< 0.89 



Taniguchi 1992; Prestwich, Joseph, & Wright 1994; Armus 
et al. 1995; Larkin et al. 1995, Scoville et al. 1998), com- 
bined with the inferred presence of ionizing sources there 
(Larkin et al. 1995; Kim ct al. 1995; Goldader et al. 1995), 
indicate that the continuous models offer a more plausible 
description of the star formation in this well studied ULIRG; 

(ii) the stellar initial mass function (IMF). The precise 
form of the IMF in starburst galaxies is very much an open 
question: see Leitherer (1998) for a review. It appears that 
the IMF of the most luminous star clusters in the Milky Way 
and the Magellanic Clouds follows closely the Salpeter form 
(Hunter et al. 1995; Hunter et al. 1997; Massey et al. 1995a). 
However, the field- star IMF may be significantly steeper 
than the cluster IMF (Massey et al. 1995b), and there is 
some evidence (Meurer et al. 1995) that the most massive 
stars in starbursts form in environments that are more sim- 
ilar to the field than to the cores of the luminous star clus- 
ters studied above. There are hence few direct observations 
to guide us to the correct IMF in starburst galaxies, and we 
therefore consider four different scenarios in Models III to 
X. Leitherer & Heckman (1995) consider power-law IMFs, 
and we investigate various possible combinations of the up- 
per and lower mass cutoffs and the slope of this power law. 
Note that if the SCUBA sources end up as elliptical galaxies 
(Section 3.1.3; Blain et al. 1999a; Lilly et al. 1998), a lower 
limit to the IMF at M\ = 3M© is suggested, in agreement 
with the value suggested by Zepf & Silk (1996). 

(iii) the inital gas metallicity. We assume a metallicity 
of twice solar. This is close to the metallicities of the most 
luminous elliptical galaxies (e.g. Faber 1973, Vader 1986); 

(iv) the age at which we observe the star clusters. We 
investigate the cases where we observe the galaxies at be- 
tween 10 8 and 10 8 ' 5 years after the onset of star formation. 
Fortunately, most of the properties of the starburst, like the 
total bolometric luminosity, reach a plateau soon after the 
onset of star formation, and so our results do not depend 
strongly on this parameter. The median total amount of gas 
consumed by a L* galaxy 10 8 ' 5 years after the onset of star- 
formation is 7 x 10 10 Mq in our models. This is consistent 



with the progenitor galaxies to the SCUBA sources being 
gas-rich normal galaxies; 

(v) the effects of extinction. Since we are only attempting 
to model the far-infrared and submillimetre properties of the 
starburst, we can assume that the galaxies are optically thin 
and do not need to make any corrections for extinction. We 
assume that the dust absorbs and reradiates all the energy 
produced by the starburst. 

From the seventh column of Table 5 it is immediately 
apparent that there is more than an order of magnitude 
of uncertainty in the transformation from far-infrared lumi- 
nosity to star-formation rate, due to the uncertainty in the 
star-formation parameters. The constant conversion factor 
of 2.2 x 10 9 L M" 1 yr (Rowan- Robinson et al. 1997), used 
by Blain et al. (1999a), corresponding to 0.22 in the seventh 
column of Table 5, is bracketed by our results. 

3.1.2 The Madau plot 

We combine our redshift-dependent luminosity functions 
(Models B, F and G) with the star-formation parameters for 
the models listed in Table 5 (Models I to X) to compute the 
comoving star formation rate in the SCUBA sources. The 
results are listed in Table 6. We can then put the SCUBA 
sources on the the Madau plot: see Figs 9, 10 and 11. The 
other points on the Madau plot all come from optically- 
selected samples. As there is no luminosity evolution in the 
models between z\ and 22, the comoving star-formation rate 
is constant between these redshifts, and so the lines on the 
Madau plot are horizontal. 

The figures show that the uncertainty in determining 
where the SCUBA sources lie on the Madau plot is huge, 
more than an order of magnitude, even if we have specified 
the correct luminosity-redshift distribution model B, F or 
G. This uncertainty is solely due to our lack of knowledge re- 
garding the star-formation parameters, as discussed in Sec- 
tion 3.1.1. For Model B, and possibly Model F, if the true 
cosmic star-formation rate in the SCUBA sources is towards 
the higher end of our permitted range, then they dominate 
the star-formation rate of the Universe: this is essentially 
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Table 4. Predicted 850-/im background intensities. The observed 
850-/im background intensity is 0.55 ± 0.15 nW m~ 2 sr _1 (Fixscn 
et al. 1998). 



120 



Figure 8. The lower limit on the integral in equation (5) as a 
function of the galaxy dust temperature T, normalized to its value 
for T = 38 K, predicted from IRAS and ISO counts by Blain et 
al. (1999a). 



the scenario proposed by Blain et al. (1999a) and Hughes 
et al. (1998). If the true cosmic star- formation rate in the 
SCUBA sources is towards the lower end of our range, then 
they do not contribute significantly to the star-formation 
rate of the Universe at any redshift. In general, the instan- 
taneous models predict higher star-formation rates. This is 
because more stars need to be formed to produce a given 
bolometric luminosity we observe more than 10 8 years after 
the starburst. 

It is interesting to note that the IMF proposed by Zepf 
& Silk (1996) for elliptical galaxies, if valid for the SCUBA 
sources, results in them being at the extreme low end of 
our proposed range. This is because these models underpro- 
duce low-mass stars, and so result in a low star-formation 
rate for a given amount of bolometric luminosity: the low- 
mass stars, when young, do not contribute significantly to 
the bolometric luminosity. The implication here is that if 
SCUBA sources evolve into elliptical galaxies, then they do 
not contribute significantly to the star-formation rate of the 
Universe at any redshift given our models of the luminosity 
function. 

The broadband optical colours of the SCUBA sources 
are not significantly different from those of normal field 
galaxies (Smail et al. 1998). Furthermore, two of the three 
local ULIRGs studied at ultraviolet wavelengths with HST 
by Trentham, Kormendy & Sanders (1998) would also have 
broadband colours similar to normal field galaxes if they 
were placed at high redshift. Therefore, SCUBA sources can- 
not be identified as submillimetre-luminous galaxies based 
on broadband optical colours alone. The star-formation rates 
given in Table 5 are far higher than the rates in normal 
galaxies. Therefore we need to treat the SCUBA sources as a 



Model 


Background / 


Model 


Background / 




nW m~ 2 sr~ 1 




nW m~ 2 sr _1 


A 


3.40 


E 


1.71 


B 


0.85 


F 


0.71 


C 


1.32 


G 


0.34 


D 


1.19 







separate population on the Madau plot. This is particularly 
important if the SCUBA sources contribute significantly to 
the star-formation rate of the Universe at any redshift. 



3.1.3 Why do the result differ from those of other authors 

Both Blain et al. (1999a) and Hughes et al. (1998) derive 
global star-formation rates at least an order of magnitude 
higher than we derive in the current paper. There are two 
sources of discrepancy. 

(i) differences in star-formation rates. For example Blain et 
al. (1999a) use a value of M\ below 1 M and consequently 
assume star-formation rates that are towards the upper end 
of our considered range. 

(ii) differences in the assumed luminosity- redshift distribu- 
tion, which can be broadly split further into (a) differences 
in the characteristic luminosity of the sources who dominate 
the contribution to the SCUBA counts , and (b) difference 
in the normalization of the luminosity function at a given 
redshift. Effect (iia) is not a major source of discrepancy 
between our results and those of Blain et al. (1999a) — in 
both cases, most of the star formation occurs in very lu- 
minous sources with L 6 o > 10 11 L . Effect (iib) is a much 
more substantial source of discrepancy. Blain et al. (1999a) 
assume a redshift-dependent Saunders luminosity function 
at high redshift which is related to the local 60- /xm lumi- 
nosity function by simple luminosity evolution. We assume 
a Gaussian luminosity function with no evolution over some 
specified redshift range. Our objects, when evolved to, z — 
are early-type stellar populations, quite unrelated to the ob- 
jects which dominate the local 60-/im luminosity function. 
We consequently derive a far lower normalization to our lu- 
minosity functions than do Blain et al. (1999a). The present 
observations do not permit us to distinguish between these 
models, but this should be a straightforward exercise when 
a significantly complete redshift distribution for the SCUBA 
sources is known. 

In addition, Eales et al. (1998) argue that at least 10 
percent of the stars in the Universe formed in SCUBA 
sources since they produce about 10 percent of the extra- 
galactic background light. This number is much closer to the 
numbers presented in this paper. Nevertheless, the assump- 
tion that the SCUBA sources contribute the same fraction of 
the submillimeter background at all wavelengths, on which 
the Eales et al. calculation is based, appears to be inconsis- 
tent with the formulation that we use, when the BIMA 2.8- 
mm counts, ISO 175-pm counts, and COBE 850-/im back- 
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Table 5. Properties of star-formation models, a is the slope of the stellar IMF and M u and M\ are the lower and upper mass cutoffs. The 
final two columns are derived from the results of Leithercr & Heckman (1995). The star-formation rates (SFR) as a function of 60-/im 
luminosity come from Figs 7 and 8, assuming a temperature of 70 K. The metal masses Mz come from Figs 53 and 54 of Leithercr & 
Heckman (1995), assuming that the mass of metals produced is equal to the mass returned by winds and supernovac. All models assume 
an initial metallicity of twice solar. 



Model 


SF profile 


logi (agc/yr) 


a 


M u / M 


Mi / Mq 


SFR / (L 60 /10 9 L Q ) M yr-! 


M z 1 (L 6 o/1O 9 L ) M 


SF I 


Instantaneous 


8.0 


2.35 


100 


1 


1.0 


4.0 x 10 7 


SF II 


Instantaneous 


8.5 


2.35 


100 


1 


0.81 


1.0 x 10 8 


SF III 


Continuous 


8.0 


2.35 


100 


1 


0.11 


3.5 x 10 6 


SF IV 


Continuous 


8.5 


2.35 


100 


1 


0.072 


7.2 x 10 6 


SF V 


Continuous 


8.0 


2.35 


30 


1 


0.15 


4.1 x 10 6 


SF VI 


Continuous 


8.5 


2.35 


30 


1 


0.12 


1.0 x 10 7 


SF VII 


Continuous 


8.0 


3.3 


100 


1 


0.34 


8.0 x 10 5 


SF VIII 


Continuous 


8.5 


3.3 


100 


1 


0.25 


1.6 x 10 6 


SF IX 


Continuous 


8.0 


2.35 


100 


3 


0.056 


3.5 x 10 6 


SF X 


Continuous 


8.5 


2.35 


100 


3 


0.043 


7.2 x 10 6 




Figure 9. The comoving star-formation density of the Universe 
contained in the SCUBA sources when their luminosity-redshift 
distribution is as in Model B. The eight dashed lines represent the 
SFR histories for models (in ascending order) X, IX, IV, III, VI, 
V, VIII, VII, II and I. The other points come from: filled triangle 

- Gallcgo ct al. (1996); open triangle - Treyer et al. (1998); open 
circle - Tresse & Maddox (1998); stars - Lilly ct al. (1996); open 
hexagons — Hammer & Florcs (1998); filled squares - Connolly 
ct al. (1997); filled circles - Madau et al. (1996); open squares 

- Pcttini et al. (1998; these include a global correction for dust 
extinction). The recent work of Glazcbrook et al. (1999) suggests 
a SFR almost coincident with the z = 0.85 point of Hammer & 
Flores (1998). 

ground are all considered in conjuction. Furthermore, it is 
unclear how we should relate star-formation rates derived 
from optical or ultraviolet luminosities of normal field galax- 
ies to star-formation rates derived from submillimetre fluxes 
of the SCUBA sources. 



Figure 10. As Figure 9, but for model F. The symbols have the 
same meanings and the lines are in the same order. 



3.1.4 The fate of the star-forming galaxies 

We have hinted that the SCUBA sources could evolve into 
elliptical galaxies. Local ULIRGs have gas densities that 
are similar to the core stellar densities in elliptical galax- 
ies (Kormendy & Sanders 1992, Doyon et al. 1994). If the 
SCUBA sources have similar morphologies to local ULIRGs, 
we might then expect them to evolve into elliptical galaxies. 

In Table 7 we present the density parameter in stars 
produced by the SCUBA sources, in our luminosity-function 
and star-formation models. For most models, particularly 
those with IMFs appropriate to elliptical galaxies (mod- 
els IX and X), these numbers are very small as compared 
with the stellar density contained in the local spheroid stel- 
lar population fl Bp h, that is elliptical galaxies and bulges: 
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Figure 11. As Figure 9, but for model G. The symbols have the 
same meanings and the lines are in the same order. 



n sph = 0.0036 (Fukugita, Hogan & Peebles 1998). In Models 
IX and X, between 1 and 4 per cent of fi sp h formed in the 
luminous high-redshift SCUBA sources, and slightly more 
if a solar-neighbourhood IMF (Models III and IV) is as- 
sumed. These numbers are sufficiently low that a scenario 
in which all the SCUBA sources evolve into elliptical galax- 
ies is consistent with galaxy formation models that predict 
low-redshift elliptical formation (Kauffmann 1996; Kauff- 
mann, Chariot & White 1996; Kauffmann & Chariot 1998). 
They are also low enough to be consistent with the observed 
paucity of red galaxies in the Hubble Deep Field and the 
consequent interpretation that less than 30 per cent of field 
ellipticals formed at high redshift (Zepf 1997; Barger et al. 
1998b). 

Various lines of argument suggest that the elliptical 
galaxies in clusters are very old. Not only did the stars form a 
long time ago (Ellis et al. 1997; Stanford et al. 1998; Kodama 
ct al. 1998), but it appears that in the richest clusters the 
most luminous elliptical galaxies themselves were assembled 
by z ~ 1 (Trentham & Mobasher 1998). It is further possible 
that most 3C radio galaxies are cluster ellipticals in the pro- 
cess of formation (Best, Longair & Rottgering 1998). This 
leads to a natural question: are the SCUBA galaxies forming 
cluster ellipticals? Low-redshift cluster ellipticals are highly 
clustered, with a bias of about 4, and an even greater bias 
at high redshifts (Fry 1996; Mo & White 1996). In compari- 
son, the reasonably uniform detection rate of SCUBA galax- 
ies in the fields studied by Smail et al. (1998) suggests that 
this is unlikely. These fields contained low-redshift clusters 
(which magnify the background SCUBA sources through 
gravitational lensing), but these low-redshift clusters are un- 
related to the hypothesized ellipticals in formation that we 
are discussing here. Instead, we propose that the SCUBA 
sources are a forming trace population of field ellipticals. 
This is consistent with the constraints outlined in the pre- 
vious paragraph. This is not to say that if we happened to 



point SCUBA at a cluster in formation we would not detect 
a large number of galaxies. 

It has been recently recognized that lyman-break galax- 
ies are highly clustered at high redshift (Pettini et al. 1998), 
and so it is worth investigating a possible link between these 
objects and the SCUBA sources. Looking at the spectra 
of the two kinds of galaxies, we see that that lyman-break 
galaxies and the SCUBA sources are very likely two com- 
pletely separate populations of objects (for example com- 
pare the 0.5 jj,m spectrum of the lyman-break galaxy 1512- 
cB58 at z = 2.7 taken by Pettini et al. 1998 using LRIS 
on Keck II with the 0.5 /tm spectrum of the SCUBA source 
SMM J02399— 0136 at z = 2.8 taken by Ivison et al. 1998 us- 
ing the MOS spectrograph on the CFHT). However, while 
this observation is interesting, it does not rule out a hy- 
pothetical link since it is possible that one type of galaxy 
evolves into the other on short timescales. Measuring cor- 
relation functions would be the best way of testing such 
a hypothetical link. Unfortunately, the forty or so SCUBA 
sources found so far do not have redshift measurements, so 
such an analysis is not possible at present. Such a calcula- 
tion will, however, be possible once a large sample of SCUBA 
sources with measured redshifts is available. 

If the elliptical galaxies that went through a 
submillimetre-luminous phase represent a sub-population, 
then their star-formation history must be similar to that of 
the rest of the remaining field elliptical population. Other- 
wise a huge scatter in the stellar mass-to-light versus light 
correlation (Oegerle & Hoessel 1991) would be produced. 

A few ellipticals (e.g. Stockton et al. 1995, Dunlop et 
al. 1996) must have formed at high redshift and do not lie in 
very rich X-ray clusters. It is not obvious that these galax- 
ies will evolve into cluster ellipticals. The existence of such 
objects lends additional support to a scenario in which the 
SCUBA sources evolve into a trace sub-population of field 
elliptical galaxies. Furthermore, many high-redshift quasars 
are surrounded by heavily enriched gas, and so there must 
have been considerable star formation in their host galaxies 
at high redshift (Matteucci & Padovani 1993). 

3.1.5 Metal production 

In Table 8, we present the total amount of metals produced 
by the SCUBA sources for our models. The average metal- 
licity of the Universe at z = is poorly known, but is proba- 
bly about half solar - see for example Section 5 of Fukugita 
et al. (1998) - if the solar metallicity is 0.0189 (Anders & 
Grevesse 1989), and is much larger than the values listed 
in Table 8. Hence, it appears that a sufficiently small frac- 
tion of cosmic star-formation activity occurs in the SCUBA 
sources that they still may be a sub-population of elliptical 
galaxies undergoing a monolithic collapse at high redshift, 
without significantly overproducing the metal mass density 
of the Universe. On the other hand, if all elliptical galaxies 
formed via monolithic collapse at high redshift, then this 
would result in a metal mass density of the Universe that 
is higher than that observed - the "cluster-field analogy" of 
Madau 1997: see also Mushotzky & Loewenstein (1997) and 
Renzini (1997). 

For comparison, if the enrichment of the Lya forest is re- 
garded as a signature of cosmic enrichment, then the metal- 
licity of the intergalactic medium at z ~ 3 is about 10~ 2,5 
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Table 6. Comoving star-formation rates in units of Mg yr _1 Mpc~ 3 in the three well fitting models of distant ULIRGs. 



LF Model 


SF I 


SF II 


SF III 


SF IV 


SF V 


SF VI 


SF VII 


SF VIII 


SF IX 


SF X 


B 


0.097 


0.078 


0.011 


0.0069 


0.014 


0.012 


0.039 


0.025 


0.0064 


0.0041 


F 


0.042 


0.034 


0.0046 


0.0030 


0.0060 


0.0050 


0.017 


0.010 


0.0027 


0.0018 


G 


0.045 


0.037 


0.0050 


0.0032 


0.017 


0.0054 


0.018 


0.011 


0.0030 


0.0019 



Table 7. The present-day values of the density parameter in stars produced by each well fitting luminosity function (LF) model, in 
each star-formation model. 



LF Model SF I SF II SF III SF IV SF V SF VI SF VII SF VIII SF IX SF X 



B 0.0010 0.00082 0.00012 7.3 x 10~ 5 0.00015 0.00013 0.00041 0.00026 6.7 x 10~ 5 4.3 x 10~ 5 

F 0.00080 0.00064 8.7 x 10~ 5 5.7 x 10~ 5 0.00011 9.5 x 10~ 5 0.00032 0.00019 5.1 x 10~ 5 3.4 x 10~ 5 

G 0.0019 0.0016 0.00021 0.00014 0.00072 0.00023 0.00077 0.00047 0.00013 8.0 x lO" 5 



to 10~ solar (Cowie et al. 1995; Songalia & Cowie 1996; 
Songaila 1997; Cowie & Songaila 1998), and may be higher 
at lower redshift. This value is close to the metallicities in 
Table 8. But this similarity is perhaps not too meaningful be- 
cause one would expect that most of the metals produced in 
the SCUBA galaxies would remain bound within the galax- 
ies. This is because ellipticals have deep gravitational poten- 
tial wells and thus low ejection efficiencies (Nath & Chiba 
1995); based on their gas consumption alone, the SCUBA 
sources must have masses in excess of 10 1 Mq. Hence, only 
a small fraction of the metals should leak out into the in- 
tergalactic medium, and so the values listed in Table 8 arc 
best interpreted as upper limits when comparing with the 
intergalactic medium metallicity, if the SCUBA sources are 
all at redshifts greater than three. We also note that most 
low-redshift elliptical galaxies have metal-poor interstellar 
media: this is not what we would expect if the metals pro- 
duced in the early history of these galaxies remained bound 
to the galaxy, as we argue above. One possible explanation 
of this apparent discrepancy is that the elliptical galaxies 
accumulate substantial X-ray halos (e.g. Forman, Jones & 
Tucker 1985), which then strip the metals out of the cen- 
tres of galaxies. Another possible explanation is that the 
SCUBA sources evolve into the few local ellipticals that do 
have metal-rich interstellar media, for example Centaurus 
A. 



3.2 The SCUBA sources as AGN 

An alternative interpretation is that the SCUBA sources 
derive their bolometric luminosities from dust-enshrouded 
AGN and not from star formation at all. A number of lines of 
argument suggest that this may well be a plausible scenario. 

First, the one SCUBA-selected source that has been 
studied in detail (SMM J02399— 0136; Ivison et al. 1998) 
contains a powerful dust-enshrouded AGN. Note however 
that it is still possible for most of its bolometric luminosity 
to be produced in a starburst. The hyperluminous gravi- 



tationally lensed galaxy IRAS F10214+4724 at z = 2.29 
(Rowan- Robinson et al. 1991; Lacy et al. 1998 and refer- 
ences therein) also contains a dust-enshrouded AGN. The 
other two known ultraluminous high-redshift dusty galaxies 
H 1413+117 (Barvainis et al. 1995) and APM 08279+5255 
(Irwin et al. 1998; Lewis et al. 1998) also contain power- 
ful AGN. If these objects are at all typical of high-redshift 
submillimetre-luminous sources, then it may well be plausi- 
ble that a signficant fraction of their bolometric luminosities 
are powered by AGN as opposed to starbursts. 

A second line of argument is based on the observa- 
tion that in local Universe, massive dark objects (MDOs) 
- probably supermassive black holes - are inferred to re- 
side in the centers of most nearby galaxies (Magorrian et al. 
1997), based on gas- and stellar-dynamical measurements 
(Kormendy & Richstone 1995). These MDOs are plausibly 
dead quasar engines that have run out of fuel. 

Magorrian et al. (1997) show that the mass of the in- 
ferred MDO is well correlated with the luminosity of the 
spheroid stellar population in the host galaxy. The total 
mass density in MDOs in the local Universe is about 1.7 
x 10 6 M Q Mpc" 3 (Haehnelt et al. 1998), about an order of 
magnitude greater than the inferred mass density of black 
holes that is required to explain the blue light of quasars 
(Phinney 1998, Haehnelt et al. 1998), unless the accretion 
efficiency in quasars is much less than 1 per cent. 

One interesting idea suggested by Haehnelt et al. (1998) 
is that most accretion onto supermassive black holes in the 
centres of galaxies is obscured by dust. If this is the case, 
then the presence of a population of dust-obscured AGN can 
explain the above discrepancy. The SCUBA sources would 
be natural candidates for this population. We use our lu- 
minosity functions to calculate the comoving mass density 
of MDOs, Pmdo required to generate the observed bolomet- 
ric luminosity of the SCUBA sources, assuming accretion at 
some efficiency e - this is the accretion luminosity divided 
by the Eddington luminosity: 
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Table 8. Cosmic enrichment from each well fitting model. The results are in units of the solar metallicity, assumed to be 0.0189 (Anders 
& Grcvesse 1989). The redshift range is specified in the luminosity function model: see Table 2. All the enrichment due to distant ULIRGs 
is completed by the lower limit to this redshift range. 



LF Model z range SF I SF II SF III SF IV SF V SF VI SF VII SF VIII SF IX SF X 



B 3 < z < 5 0.12 0.29 0.010 0.021 0.012 0.030 0.0023 0.0047 .010 0.021 

F 2 < z < 4 0.050 0.13 0.0044 0.0090 0.0051 0.012 0.0010 0.0020 0.0044 0.0090 

G 1 < z < 3 0.055 0.14 0.0048 0.0098 0.0056 0.014 0.0011 0.0022 0.0048 0.0097 



Pmdo = 10 4 q (^j- j M MpcT 3 . (9) 

q = 3.8, 1.6 and 1.8 in Models B, F and G respectively. 
These numbers are all below the local mass density of 
MDOs. Therefore it is plausible that all the SCUBA sources 
are AGN, but they may only help to explain the discrep- 
ancy between the local density of MDOs and the inferred 
density of black holes required to explain the blue light of 
quasars if their accretion luminosity is significantly lower 
than 1 per cent of the Eddington value. Recent work by Al- 
maini, Lawrence & Boyle (submitted), based on the X-ray 
luminosity function and background confirms that a large 
fraction of the SCUBA sources could be dust-enshrouded 
AGN. 

A third, but weaker, line of reasoning is based on 
the analogy with the ultraluminous z — 0.44 cD galaxy 
IRASP09104-4109 (Kleinmann et al. 1988), which with a 
bolometric luminosity of about 3 x 10 13 L<0 is by far the most 
luminous object detected at mid-infrared wavelengths at 
z < 0.5. It is at the centre of a very substantial cooling flow 
(Fabian & Crawford 1995). Two of the cD galaxies in the 
lensing clusters studied by Smail et al. (1998) were detected 
at 850 lira (Edge et al. submitted) , although they are not in- 
cluded in the counts presented by Smail et al. (1997; 1999). 
If the physical processes that make IRAS P09104— 4109 into 
a hyperluminous galaxy are related to its being in the cen- 
ter of a cluster, then we might expect the same processes 
in the cDs detected by SCUBA. The bolometric luminos- 
ity of IRAS P09104— 4109 appears to come from a power- 
ful dust-enshrouded AGN (Kleinmann et al. 1988), and an 
analogous mechanism could plausibly be powering the high- 
redshift SCUBA sources. 



4 LIMITATIONS OF THE PRESENT WORK 
AND FUTURE PROSPECTS 

The main limitations of the present work are that: 

(i) we cannot determine whether the SCUBA sources are 
star-forming galaxies or AGN. Clearly the implications for 
cosmology would be very different for these two scenarios; 

(ii) wc cannot determine uniquely the redshift-dependent 
luminosity function of the SCUBA sources, either to tell 
which of Models B, F and G or a different kind of model 
(Blain et al. 1999a) are most correct; 

(iii) we cannot determine the conversion factor to trans- 
form the far-infrared luminosity into a star-formation rate 
for the SCUBA sources to within one order of magnitude, 



even if all the far-infrared luminosity is assumed to be gen- 
erated by star formation. 

Given these limitations, wc suggest the following future 
studies to address each of the problems listed above: 

(i) Optical spectroscopic observations of SCUBA galaxies 
should be made to search for AGN, for example to determine 
the line ratios of hydrogen recombination to forbidden lines 
(Sanders 1998b; Veilleux et al. 1995; Evans 1996) and carry 
out spectral polarimetry (e.g. Hines & Wills 1993, Januzzi 
et al. 1994). Note, however, that these measurements are ex- 
tremely difficult to carry out for very distant galaxies, and 
cannot be used to determine unambiguously whether most 
of the bolometric luminosity comes from a starburst or an 
AGN, because of the huge amount of internal extinction 
towards the galaxy centres at optical wavelengths. Genzel 
et al. (1998) used the ratio of high- to low- excitation mid- 
infrared emission lines in combination with the strength of a 
7.7-/xm PAH feature to search for AGN in local ULIRGs, but 
it will be quite some time before such measurements become 
available for the distant SCUBA sources. Crude interpreta- 
tions can be made on the basis of source temperatures and 
SEDs (Sanders et al. 1998a), for example compare the SEDs 
of Arp 220 and Mrk 231. Perhaps the best we can hope for 
in the near future is to match the SED of each SCUBA 
source to that of a plausible low- redshift counterpart, for 
which mid-infrared observations have been used to deter- 
mine a power source, and assume the high-redshift source is 
powered in the same way; 

(ii) an accurate determination of the redshift-dependent 
luminosity function will become available once redshifts 
for a significant number of SCUBA sources are ob- 
tained. In addition, further surveys with SCUBA and other 
millimetre/submillimetre-wave detectors (Glenn et al. 1998) 
should allow the luminosity function to be determined accu- 
rately (Blain 1998). For many of our models (Table 8), the 
total star-formation rate in high-redshift ULIRGs is greater 
than at low redshifts; a measurement of the redshift distri- 
bution would allow this statement to be quantified far more 
accurately. 

(iii) selecting star-formation parameters to use in the Lei- 
therer & Heckman models is a difficult exercise. This is 
a particular problem for the distant SCUBA galaxies, for 
which only very limited information about their global prop- 
erties is available. Progress in this area would be possible if 
the low-redshift remnants of the SCUBA galaxies were iden- 
tified, and their stellar populations studied in detail in order 
to investigate their IMF. Some clues are given in Section 3.1, 
but the conclusions there are far from rigorous at this stage. 
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More generally, placing the SCUBA galaxies in the con- 
text of the galaxy-formation problem is clearly an important 
exercise. Measurements of the redshift distribution will be 
an crucial part of this. For example, these measurements 
will let us see where the comoving bolometric luminosity 
density of the SCUBA sources peaks, relative to the peak in 
the comoving star-formation rate in other galaxies (Madau 
et al. 1996, Madau et al. 1998, Pettini et al. 1998) and to 
the peak in the comoving AGN luminosity density at ul- 
traviolet (Boyle & Terlevich 1998) or radio (Dunlop 1998) 
wavelengths. 



5 CONCLUSIONS 

(i) We identify various redshift-dependent luminosity 
functions that are consistent with the 850-/xm SCUBA 
source counts and observations at other far-infrared and sub- 
millimetre wavelengths. Consistency can only be achieved if 
the 850- fim SCUBA source counts are towards the lower end 
of the range suggested by Smail et al. (1997). In all of the 
consistent models, the SCUBA sources have high bolometric 
luminosities in excess of 10 12 Lq. 

(ii) At present we cannot identify whether the power 
sources in the SCUBA galaxies are predominantly starbursts 
or AGN. 

(iii) If the SCUBA sources are dusty ultraluminous star- 
forming galaxies, then they need not dominate the star- 
formation rate of the Universe at any redshift. In particular, 
if they evolve into elliptical galaxies, our models suggest that 
they only contribute a few percent of the local spheroidal 
stellar population. This result is consistent with a scenario 
in which most elliptical galaxies formed at low redshift, as is 
implied by the paucity of red galaxies in wide-field AT-band 
surveys. 

(iv) The alternative scenario in which the SCUBA sources 
are powerful dust-enshrouded AGN is also possible. If the 
AGN are powered by accretion onto a black hole with a lu- 
minosity that is about one per cent of the Eddington value, 
then they may help to explain away the discrepancy of a 
factor of ten between the observed mass density of super- 
massive black holes in the local Universe and that required 
to produce the blue light of high-redshift quasars. 



(v) r . 'ight constraints on the redshift-dependent luminos- 
ity function of submillimetre-luminous galaxies at high red- 
shift will follow from spectroscopic measurements of a large 
sample of SCUBA sources (Barger et al. in preparation). 
Obtaining accurate star-formation rates for the galaxies and 
identifying the nature of the power source will be more dif- 
ficult. Matching the SEDs of the SCUBA sources to those 
of local ULIRGs, for which power sources have been deter- 
mined using mid-infrared spectroscopy, offers one way to 
make progress in this regard. 
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